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Abstract. The Cauchy-Schlomilch transformation states that for a function 
/ and a, b > 0, the integral of f{x^) and af{{ax — bx~^)-^ over the interval 
[0, oo) are the same. This elementary result is used to evaluate many non- 
elementary definite integrals, most of which cannot be obtained by symbolic 
packages. Applications to probability distributions is also given. 



1. Introduction 

The problem of analytic evaluations of definite integrals has been of interest to 
scientists for a long time. The central question can be stated as follows: 

given a class of functions ^ and an interval [a, b] C M, express the integral of f Cz ^ 



J a 

in terms of special values of functions in an enlarged class ©. 

Many methods for the evaluation of definite integrals have been developed since 
the early stages of Integral Calculus, which resulted in a variety of ad- hoc techniques 
for producing closed-form expressions. Although a general procedure applicable to 
all integrals is undoubtedly unattainable, it is within reason to expect a systematic 
cataloguing procedure for large groups of definite integrals. To this effect, one of 
the authors has instituted a project to verify all the entries in the popular table by 
I. S. Gradshteyn and I. M. Ryzhik [13]. The website 

http : //www .math, tulane . edu/~vhm/web_ html/pap- index .html 

contains a series of papers as a treatment to the above- alluded project. 

Naturally, any document containing a large number of entries, such as the table 
[13! or the encyclopedic treatise |20j . is likely to contain errors, many of which 
arising from transcription from other tables. The earliest extensive table of integrals 
still accessible is [2] , compiled by Bierens de Haan who also presented in [3] a survey 
of the methods employed in the verification of the entries from [T3|. These tables 
form the main source for |13| . 

The revision of integral tables is nothing new. C. F. Lindman [17] compiled a 
long list of errors from the table by Bierens de Haan [4] . The editors of [T^ maintain 
the webpage 

http : //www . mathtable . com/gr/ 
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where the corrections to the table are stored. 

Many techniques have been developed for evaluating definite integrals, and the 
goal of this paper is to present one such method popularized by Schlomilch in 
|24) . The identity (12. ip below appeared in where it was called the Schlomilch 
transform, although it was used by J. Liouville [18] to evaluate the integral 



(1.1) 



1 iM+i/2(i_^y.-i/2^^ 
{a + bt~ cf^Y+^ ' 



and in ITOl Liouville quotes a letter from Schlomilch in which he describes his 
approach to (|l.ip via the formula 

to derive the reduction 



(I 3) 

Iq (a + ^o; + 7x2)"+! ^ J„ (^ + 2^57 + y)'"+i ' 

The integral ()l.ip is then evaluated in terms of the beta function. Schlomilch also 
states that his method can be found in a note by A. Cauchy |8j published some 35 
years earlieiQ. In view of this historical precedence, the name Cauchy-Schldmilch 
adopted here seems to be more justified. Some illustrative examples appear in the 
text [21]. 

We present here a variety of definite integrals evaluated by use of the Schlomilch 
transform and its extensions. Several of the examples are not computable by the 
current symbolic languages. For each evaluation presented in the upcoming sec- 
tions, we considered its computation using Mathematica. Naturally, the statement 
'the integral can not be computed symbolically' has to be complemented with the 
phrase 'at the present date'. 

2. The Cauchy-Schlomilch transformation 

In this section we present the basic result accompanied with initial examples. 
Further applications are discussed in the remaining sections. 

Theorem 2.1. [Cauchy-Schlomilch] Let a, b > and assume that / is a continuous 
function for which the integrals in ()2.1|) are convergent. Then 



/•oo -1 /'OO 

(2.1) f (^{ax - bx-^y) dx = -j^ f{y^)dy. 

Proof. The change of variables t — bj ax yields 

1=1 f (j^ax — b / x)^^ dx 



b 







f[{at^b/tf) t- 



dt. 



The average of these two representations, followed by the change of variables u = 
ax — b/x completes the proof. □ 

The next result is a direct consequence of Theorem 12.11 



^This note is available at http : //gallica.bnf . fr/ark: /12148/bpt6k90193x. zoom . f 526 
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2fe+l 

k=0 



Corollary 2.2. Preserve the assumptions of Theorem 12.11 Let /i„ ~ CkX 
be an odd polynomial. Define 

(2.2) g„(x) j^dfcxM , 

\fc=o / 

where 



j=k 

Then 

/•OO 1 /"OO 

(2.4) / ([/^„(ax) - h„{bx-'f) dx=-J^ f{9n{y')) dy. 
Proof. Denote 

n 

(2.5) Hn{x) = hn{ax) - hn{bx^^) = Cklpkjx) 

where V'fc(a;) = (aa;)^''+^ - (62;"^)^''+^ and let <j)k{x) = {ax - bx^^)'^''+'^ . Then, the 
polynomials ^k and (f>k obey the transformation rule 

(2.6) M^) = E 

The expression for iJ„ in terms of (pk follows directly from this. Moreover, H^{x) — 
gn {{ax — bx^^)'^) and the result is obtained by applying the Cauchy-Schlomilch 
formula to f{gn)- O 

Example 2.3. For n e N, 

r [(.^ + + 1, - 1] e-.-'-'.- J-r (i.) . 

In order to verify this value, write the Laurent polynomials of the integrand ac- 
cording to (|2.3I) 

- a;""^ ^7y + Uy^ + 7y^ + y"^, 
and write the integrand as a function of y — x — x^^ , to obtain 

7{x^ + x-'^){x'^ ^x^ + l)-7 = 7 + 42y2 + 35y4 + 7y^ 
Implement (|2.4p on 

(2.7) f{y^) = {7 + 42y2 + 35/ + 7y6)g-(y^(7+i4^;^+7/+,«)^)" ^ 
to yield 

poo pOO 

7 [{x"^ + x-^){x'^ - + 1) - 1] e-(^'-^-')'" dx = f{y^)dy. 

Jo 

The last step is an outcome of the substitution z ^ 7y + 14y^ + 7y^ + y'^ . Hence 
/ f{y^)dy = / (7 + 42y2 + 35/ + 7/)e-(^^+"^'+^^'+^')'"dy 
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Example 2.4. Proceeding as in the previous example, we obtain 

r [7{x^ + - x2 + 1) - 6] dx = -^r (^) 

Jo 2n \2?T./ 



The details are left to the reader. 



Example 2.5. The choice y — x — x ^, x'^ — x ^ = 3y + y^ followed hy z = 3y + y^ 
produces 

- — n [1 + - ^"')' ^''] ~ 



Jo 



= -(l+i/+I/''^ + !^'^ + i^^°H 

6 

3. An INTEGRAL DUE TO LAPLACE 

The example described in this section is the original problem to which the 
Cauchy-Schlomilch transformation was applied. 

Example 3.1. The normal integralis 

(3.1) fe-y'dy^^. 

The reader will find in [6] a variety of proofs of this fundamental identity. 
Take f{x) — e^^ in Theorem 12. ll to obtain 



(3.2) 



} / — 

2a 



Expanding the integrand and replacing the parameters a and b by their square roots 
produces entry 3.325 in [13]: 

(3.3) / exp (— ax^ — 6/2:^) = ^1 /— e 

Jq 2 \ a 

The change of variable x = \/bt/^/a shows that the result p.3p can be written 
in terms of a single parameter c = a6 as 

Example 3.2. A host of other entries in [13] are amenable to the Cauchy-Schlomilch 
transformation. For example, 3.324.2 states that 



(3.5) 



This is now evaluated by choosing /(x) — e in Theorem l2.1l so that 



/oo pc 
cxp[-(.T-&/a;)^"] dx^2 
-00 Jo 



dy. 
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The change of variables t = y^" and the integral representation for the gamma 
function 



POO 



PC 

(3.6) r(a) = / e~*r'Ut 

Jo 

imply (3.5). 

Example 3.3. The expression t — 1 /t in (|3.4p suggests a natural change of variables 
i = e". This yields 

(3.7) / e"^" p. 



The latest version of Mathematica is unable to produce this result when c is an 
arbitrary parameter. It does evaluate p.7p if c is assigned a specific real value. 

4. An integral with three parameters 

The introduction of parameters in a definite integral provides a greater flexibility 
in its evaluation. Many classical integrals are presented in [5] as special cases of 
the next theorem, which now appears as 3.242.2 in [T3]. The proof given below is 
in the spirit of the original observation of A. Cauchy. 



Theorem 4.1. Let 

h 



I: 



2 



a;4^ 


h 2ax^ ^ 


h 1 










X 






1- 2aa;2 - 


h 1 










X 





1 dx 
1 ^ 



dx 
a;2 



^ Jo W + 2ax^ + l ) 

^ ^ 2Jo \x^ + 2ax^ + l) x^ 
Then Ii = I2 = k ~ li and this common value is 
(4.1) I{a, b; c) = 2-1/2-^(1 + 0)1/2-^5 - i, i) . 

Proof. Observe that if g satisfies g{l/x) — x'^g{x), differentiation with respect to 
the parameter b shows that the integral of g{x)/{x^ + 1) over [0, 00) is independent 
of b. This proves the equivalence of the four stated integrals. 

Theorem 12.11 is now used to evaluate I3. For any function /, the Cauchy- 
Schlomilch transformation gives 







„2 



^'-'2a..2 + lj^- = i H(----^)^ + 2a + 2j'^" 

/ 



1 

1 

1 



+ 2(a+ 1) 



Apply this to f{x) — x'^ and use the change of variables u — j.2'^2(a+i) produce 



6 T. AMDEBERHAN, M. L. GLASSER, M. C. JONES, V. MOLL, R. POSEY, AND D. VARELA 

This last integral is the special value B{c — ^, ^) of Euler's beta function. 

□ 

The next theorem presents an alternative form of the integral in Theorem 12. II 
Theorem 4.2. For any function /, 

where a, = 2(TT^- 

Proof. This follows from the identity in Theorem 12.11 and the change of variable 
t^2{a + l)/[x^ + 2{a + l)]. □ 

The master formula (j4.ip yields many other evaluations of definite integrals; see 
[7] for some of them. The next theorem provides a new class of integrals that are 
derived from (|4.ip . 

Theorem 4.3. Suppose 

oo 
n=l 

be an analytic function with /(O) — 0. Then 
where u = -^jy-. — r. 

8(l+a) 

Proof. Integrate term-by-term and use the value 




to simplify the result. □ 



5. Exponentials and Bessel functions 



This section describes the application of Theorem 14.31 to a number of definite 
integrals. The Taylor expansion oi f{x) — 1 — e^^^ with 5 > leads to an integral 
that can be evaluated in terms of the modified Bessel functions I^{x) defined by 
the series 

In particular 

oo 2j oo 2j + l 

(5.2) Io{x) ^ -^i(^) ^ E 22j+i^,l(, + IV • 

Although, as will be seen below, our results below have more direct derivations, the 
following procedure is more informative. 
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Example 5.1. For a > —1 and 6 > 0, let c = 8{i+a) • Then 

f (l - e-^^V(.^+2a.^+i)) ^ M2c) + hi2c)] . 

The function f{x) = 1 - e"''^ has coefficients c„ = and (gSJ yields 

V / 23/2 ^/TT^ ^ ^ 



where w = 1/8(1 + a) and 
(5.3) hix) = 



23/2 yrrp-a 

('") 



(n+l)! 

The result now follows from the relation c ^ bu and an identification of the series 
h in terms of Bessel functions. 

Proposition 5.2. The following identity holds: 

oo /2n\ 

We present two different proofs. The first one is elementary and is based on 
the WZ- method described in [53]. Mathematica actually provides a third proof by 
direct evaluation of the series. 



Proof. The expansion ()5.2p yields 

(5.4) Io{2x)-hi2x)^J2 
where 



OO ^ 

X 



r—0 



ifr = 2j 
+ ifr-2j + l. 

Multiplying (|5.4p with the series for e^^ lends itself to an equivalnet formulation of 
the claim as the identity 

^' (-ly fk\ f j \ 1 /2/^^ 



(5.5) J2 



The upper index of the sum is extended to infinity and (|5.5p is written as 

{-ly fk\ f j \ _ 1 f2k^ 



2j + l 
j 



^0 Vj7VLj72J; 2Hk + i)\k 

The even and odd indices are considered separately. Define 

The result is obtained from the values 
(5.6) ^e^iP'Und^o ' 



2'^Vfcy (A: + l)2feVfc 
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To establish (|5.6p . the WZ-method is applied to the functions 

S:{k) = 5e and S:{k) = 5. 

The output is that both S* and 5** satisfy the recurrence Uk+i — = with 
certificates 

(2fc+l)(/c + l + 2j) (2A: + l)(/c-2j)' 
respectively. The initial conditions S*{0) = S'*(0) = 1 give S*{k) = S'*(fc) = 1. □ 

As promised above we present an alternative proof of I5.2l based on Theorem 14. 2 1 
Proof. Theorem 14.21 gives 

Vb /"M-e-"** dt 



OO , 

2 //^4 



1 _ e-''-V(-*+2a.^+i) rf^^ 



^ ' Jo t VtU~~t)' 

The latter is known as Frullani integral and can be written as 

5.7 I = / dydt. 

^ ' 2 Jo Jo VtU^ 

Exchanging the order of integration, the inner integral is a well-known Laplace 
transform [11 (p. 366, 19.5.11 with n = 0) 

whence we find 

1=-^ r '\-^i,(t)dt. 

Va* Jo 

The relation 

(5.9) |(te-*(/o(t)+/i(0))=e-%W 

now completes the proof. □ 



Example 5.3. Choosing a — and b — 4 gives 

V J e 

Example 5.4. The values a — I and 6 = 8 yield 

2n 



Mathematica is unable to evaluate these two examples. 
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6. Trigonometric and Bessel functions 

The next example employs the familar Taylor expansion of sin&a;. The result is 
expressed in terms of the Bessel function of the first kind 

Example 6.1. Let c = 6/8(1 + a). Then 

f°° ( hx^ \ nb 

/ sin — ; ^ dx — — , [ Jn (2c) cos 2c + Ji (2c) sin 2cl . 

Jo \x^ + 2ax^ + lj ^8(1 + a) ^ ^ ' n 7 j 

To verify this, apply Theorem 14.31 to sin 6a; to obtain 

X' 







4 + 2ax^ + 1 ; ^8(1 + a) (2A; + 1)! V2fc 



Lemma 6.2. The following identity holds: 

\k 



JljkTYy\2k)''' = ^o(2c)eos2c+ Ji(2c)sin2c. 



Proof. Using the Cauchy product and the series expression (|6.ip . the right-hand 
side is written as 

Jo (2c) cos 2c + Ji (2c) sin 2c = 



oo k 



{-Ifc/H^ - - (-l)*c^MJ 



y-y- V-r; c I 2c^ V V 



(-l)^c^^(4fc)! - (^l)fcc^fc(4fc + 3)! 

(2^)!' ' ^t^„(2fc + 3)!(2fc+l)2 



(-I)'' /4fc 



..0 (2^ + l)!V2fc 



c2^ 



The passage from the first to the second equality is justified by the identities 

A 4^ _ 1 /4fc\ A _ 4fc + 3 /4fc + 2 

^(2j)!(fc-j)!2 " (2fc)!Uy '''' (2j)!(fc - j)! " J + 1)! " (2fc + 3)! \2k + 1 

Both of these formulas are in turn verifiable via the WZ-method t23j with their 
respective rational certificates 

(6n2 + 10n + 4-4nfc- 34A:)(2A:- l)/fc (20n + 17 + Gn^ - 4n/c - 7fc)(2fc - 1) 

RIIQ 



{n + l- /fc)2(4n + l)(4n + 3) (n + 1 - k)(n + 2 - k){An + 5)(4n + 7) ' 

□ 



Example 6.3. The choice a — Q and 6 = 1 in Example 16.11 produces 
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Example 6.4. By choosing a = 6 = 1 in Example 16.11 we get 

2 

/ / T" 

(6.3) 



a;2 + l 



<i^^^[JoCs) cos(i)+Ji(i) sin(i)] 



As the time of this writing, Mathematica is unable to evaluate the integrals in the 
two previous examples. 



Note. The function 



(6.4) 



5(") = E 



k=0 



(2fc + 1)! \2k 



,2k 



also has a hypergeometric form as 
(6.5) g{u) = 2F3 

This follows directly from the identity 



1 3 
4 4 



MfeN 
V2fe/ 



22fc-3/2r(fc + i/4)r(fc + 3/4) 



(2A: + 1)! r(/c + l/2)r2(/c + l)r(fc + 3/2)' 
that is established via the duplication formula of the gamma function 



1 



r(2x) = -^2^^-^T{x)r{x + ^) 



and its iteration 



r(4a;) - J—2^^-^/^T{x)r{x + ^)T{x + i)r(x + f ) 
ttVtt 



Then r(a + k) = {a)kT{a) produces l|6.5p . Here {a)k = a{a + 1) • • • (a + fc — 1) is 
the Pochhammer symbol. This gives an alternative form of the result described in 
Example 16. li i.e., 



bx' 



dx 



TT b 



2F3 



^a;4 + 2ax^ + 1 
Note. Proceeding as in Example 12. 31 we can obtain 
(x^ + a;^2 _ 1) gin 



1 3 
4 4 

1 1 3 

2 2 



-b^/16{l + af 



jl2 _ 42-6 _ 42,-6 



12 



dx 



6V2 



[Jo(i)cosi + Ji(i)sini] 



Mathematica is unable to compute the preceding integral. 

Note. The result in Example 16.11 can also be established by the Frullani method 
described in Section [5j Start with 



/ := 
As before, write 



sm 







bx^ 



x^ + 2ax^ + 1 



dx = 



/"^ sin(a*i) dt 



2Ja* 



I = 



cos{a*ty) 
Jo ^/til - t) 



dt dy. 
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The inner integral is a well-known cosine transform (^i)j,p.l2, 1.3. 1^1 

, ^, f^cos{ujt)dt ^ /w\ uj 



and it follows that 



'a Jo 

The result is now immediate from the identity 



cost Jo (i) dt. 



^ [t {costJo{t) + sint Ji(i))] — costJo{t). 

7. The sine integral and Bessel functions 
The sine integral is defined by 

(7.1) Siix) -.^ £ dt. 

The Cauchy-Schlomilch transformation can be employed to prove 
f°° ( hx^ \ 

/ ~ ^ — 2 7 '^^ ^ 77^2(1 + a) [(4c cos 2c- sin 2c) Jo (2c) + 4c sin 2c Ji (2c)] 

Jo V*^ ^ l(xx + ly 

where c = 6/8(l+a). To establish this identity, start with the evaluation in Example 

L "" U^ + 2ax^ + l j '^-=^|^['^o(2c)cos2c+J,(2c)sin2c] 
divide by h and integrate both sides. Then the identity 

(7.2) — \(2x cos X — sin x) Jo (a:) + 2x sino; Ji(a;)] — Jo(x) cos x + Ji(x) sin a; 
dx 

gives the result. 

Example 7.1. The special case a = and 6=1 yields 

I rf^ = ^[^o(i)[cosi-2sini]+Ji(i)sini]. 

Example 7.2. The special case a = b = 1 yields 



2\ 

I T 

Si 







x2 + 1 

Mathematica is unable to evaluate these integrals. 



^^=2 ['^°'^^) [cos|-4sini] + Ji(|)sini] 



^Note that formula 1.3.14 in |10| is incorrect. 
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8. The Riemann zeta function 

Interesting examples of definite integrals come from integral representations of 
special functions. For the Riemann zeta function 

(8-1) C(^) = E;;7 

one such expression is given by 

(8.2) C(^)= I t dt 



k=l 



il-2^-^)T{s) Jo 1 + e* 

Analytic properties of C(s) are often established via such integral formulas. 
The change of variables t — produces 

poo 2s-l J 1 

(8.3) ^^-^^-(l-2i-)r(.)C(s). 
Example 8.1. We now employ Theorem l2.1l to establish 

foo 2s+l J 

(8.4) / — 2-|- = 2-(l-2i-)r(,s + l)C(s). 
Jo cosh (X'') 

The notation 

(8.5) A(s) := i^|^r(.)C(s) 

is employed in the proof. First introduce the change of variable x = f in the 
Cauchy- Schlomilch formula and take a = h for simplicity. Then 

/ f-^f[a\f -t-^f]dt^- f{y^)dy. 
Jo ^'^ Jo 

Now let f{x) = x^-^/^/{l + e^). Using the notation Sr = sinh(rlni), (|8^ yields 

A{s) _ r°° r-^{at~ at-y'^-Ut _ /-^ r-\2aSrf'^Ut 
"oT^X 1 + exp[(ar - at-^)^] ~ 1 + exp[(2aS'r)2] ' 

Differentiate (|8.7p with respect to a and use 
c 1 



(8.6) 



(l + c)2 4cosh^(6'/2) 

to obtain 

2s-l A(^ _ 2f2a)2-i /"^ t^-'S^/+Ut ^ _ A(£) 
a ar Jo cosh2[2(a5'^)2] a^r 

that produces 

e-\S^'+^dt 2sA(s) 



(8.10) 



cosh2[2(a5,)2] (2a)2«a2^- 

Now change t by in (|8.10l) and average the resulting integral with itself. The 
outcome is written as 

r-i52-+irft _ 2gA(.) 
^ ■ ^ 7o cosh2[2(a5,)2] (2a)2«a2^- 

The final change of variables x — \/2aSr produces the stated result. 
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The special case s — ^ yields 
(8.12) 



/o cosh2(a;2) 4 
Mathematica is unable to produce (|8.12p . 



"^'^"^ =-i(2-^/2)C(l/2)V^. 



9. The error function 
Several entries in the table |13| involve the error function 

erf(a;) :— —= 

V7I- 

For instance, entry 3.466.1 states that 



(9.1) 



e"* dt. 



(9.2) 



(l-erf(^/3))e^ ^ 



x2 + ^2 2/3 

Mathematica is able to compute this example, which can be checked by writing the 
exponential in the integrand as e'' '^^ x "'"'^^ and differentiating with respect 

to 

Example 9.1. The Cauchy-Schlomilch transformation is now applied to the func- 
tion 



x"^ + 2a + x~2 
The choice a = /i = 1 yields 



(9.3) 

to produce 
(9.4) 



/(x2 + 2(a+l)) 



7re 



2^2(a+ 1) 



1-erf 



(Ai\/2(a+ 1) 



(9.5) 

and a = 0, /i = 1 gives 
(9.6) 



(x + .-^)2 = - [1 --'^f(^)^ 

" ^l-crf(y2) 



2V2 

Neither of these special cases is computable by the current version of Mathematica. 

10. Elliptic integrals 



The classical elliptic integral of the first kind is defined by 



(10.1) 



K(fc) 



dx 



r/2 



dip 



lo ^{l-x^){l-k^x^) Jo y/l - fc2 sin^ cp 

The table [13| contains a variety of definite integrals that can be evaluated in terms 
of K(A;). For instance, entry 3.843.4 states that 

tanx 



(10.2) 



= ^ ^K(fc). 
1 - fc2sin2(2a;) ^ 



The reader will find in [IB] a large variety of examples. 

In the context of the Cauchy-Schlomilch transformation, we present two illus- 
trative examples. 
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Example 10.1. The first result is 

,10.) /" ^ ^=^K 



^/{x^ + 2ax2 + l)(x4 + 2bx^ + 1) ^^2(0+1) yVa+1 
To verify this result, apply (|2.ip to the function 



^y{x + 2a + 2){x + 2b + 2) 
and observe that 



^(a;4 + 2aa;2 + l){x* + 2bx^ + 1) 
The Cauchy-Schlomilch transformation gives 

a;2 dx f°° dx 



v/(2;^ + 2ax2 + l)(a;4 + 2bx^ + 1) 7o V(a^^42)(a?2 + S2) 

with = 2(a+ 1) and = 2{b+ 1). This last integral is computed by the change 
of variable x = A tan if and the trigonometric form of the elliptic integral yields 

Example 10.2. A similar argument produces the second elliptic integral evalua- 
tion. This time it involves 

(10.4) F{^,k) ■ '^^ 



Jo y/{l-x^)il-Px^) 

the (incomplete) elliptic integral of the first kind. 
Assume a < b < c. Then 

(10.5) ^'''^ 



a/(x4 + 2ax2 + l)(a;4 + 26a;2 + l)(x4 + 2ca;2 + 1) 
1 



2y/{b + l){c-a) 



F 



sin ^ 



c-a {b-a){c+l) 
V 5+l)(c-a) 



Following [12] a generalization of the Cauchy-Schlomilch identity is now used to 
evaluate some hyper-elliptic integrals. 

Theorem 10.3. Assume (j){z) is a meromorphic function with only real simple 
poles aj with Res((/);aj) < 0. Moroever assume (f) is asymptotically linear. Then, 
for any even real valued function /, 

poo poo 

(10.6) / f[(l){x)] dx= f{x)dx. 

Jo Jo 

Example 10.4. As a simple illustration we take 

^ z2 _ ^2 

(10.7) ^^(,) = ,J|__|, 
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where ai < 02 < • • • < ajv- Take = 1 and write 61 = b. Theorem 110.31 and 



(10.8) / - - Itc ( 



give 



io l/fHijOli) a \ a J' 

with 

P(t) = t"" + (a^ - 2a^)t^ + {a^ - 2a^b^)t + aH\ 

Q(t) = + [P^ ^ 2a^)t^ + (a^ - 2l3%'^)t + . 

As an interesting special case, we have 

(10,10) /" . iK'ik), 

Jo y/t [t^ - 2{kk')H + fc2] [t3 - A{kk')H^ + /c4] fc 

where fc' is the complementary moduhis and K'(fc) K(fc'). Mathematica is 
unable to deal with this. 

11. An extension of the Cauchy-Schlomilch method 

An extension of Theorem 12.11 by Jones [Tl is discussed here. The next section 
presents statistical applications of this result. 

Theorem 11.1. Let s be a continuous decreasing function from M+ onto R+. 
Assume / is self-inverse, that is, s~^{x) = s{x) for all x G Then 

/•oo /"OO 

(11.1) / f{[x-s{x)]') dx^ / f{y^) dy, 
Jo Jo 

provided the integrals are convergent. 

Proof. The change of variables t = s{x) yields 

/■oo />oo 

(11.2) / f{[x~s{x)f)dx^- f{[s{t)~tf)s'{t)dt. 
Jo Jo 

The average of these two representations, followed by the change of variables u = 
X — s{x) gives the result. □ 

Note. The above result is given without a scaling constant a > 0. This could be 
introduced by the change of variable xi = ax in (jll.ip to obtain, after relabeling 
x\ as X, 

/•oo -1 /"OC 

(11.3) / f{[ax-s{ax)Y)dx = - f (y^) dy. 
Jo ^ Jo 

Jones [13] lists several specific forms of self-inverse s{x) along with two methods 
for generating such functions based on work of Kucerovsky, Marchand and Small 
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Example 11.2. An attractive self- inverse function is 

(11.4) s(x) = x - i log (e"^ - 1) . 

a 

Then (|ll.ip becomes 

(11.5) f (^^ log2(e"- - 1)) dx = fiy') dy. 



The choice f{x) = e ^ gives, using (|3.1 



(11.6) 



exp jlog^(e"^ -1)] dx 



Example 11.3. Several other examples of self-inverse functions are provided in 
Jones [14]. Each one produces a Cauchy-Schlomilch type integral. Some examples 
are 



f [{x - exp{a/logx)f] dx = / f{y^)dy, 



log' 



1 -I- cosh(Q;a;) 



dx 



fiy^)dy, 



/ [(x - sinh(a/sinh ^ x))^] dx ^ f{y^)dy. 
Jo 



12. Application to generating flexible probability distributions 

There has recently been renewed interest in the statistical literature in generat- 
ing flexible families of probability distributions for univariate continuous random 
variables. Baker [I] describes the use of Cauchy-Schlomilch transformation to gen- 
erate new probability density functions from old ones. Jones [14| does the same 
with the extended transformation of Section [TT] The identity (jll.3|) states that 
the total mass of /(y^) is the same as that of af{[ax — s(ax)]^) for any self-inverse 
function s and any scaling constant a > 0. 

There are many techniques for introducing one or more parameters into a sim- 
ple 'parent distribution' with probability density function g, to produce more so- 
phisticated distributions. One such method, not generally familiar, proceeds by 
'transformation of scale', defining fb{x) cx g(ti,{x)) where th{x) depends on the new 
parameter b. The difficulty associated with this procedure is the validation that /f, 
is integrable and then to explicitly provide its normalizing constant. 

The Cauchy-Schlomilch result in Theorem l2 . 1 1 guarantees that the choice tii{x) — 
\x — bx~^\ produces from the density 5 of a positive random variable a new density 
/b, also for a positive random variable, via 



(12.1) 



Mx)=g{\x-bx-'\). 



This was observed by Baker 1 . The parameter a in p.ip is redundant for distri- 
bution theory work since its action as a scale parameter is well understood; a must, 
however, be reintroduced for practical fitting of such distributions to data. 

A number of general properties of distributions with density of the form ff, follow, 
some of which are: 

(i) fb is R-symmetric [22] about R-center Vb, i.e. fbiVbx) — fb{Vbx^^)', 
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(ii) /fc(0) = and 

f},{x) w g{b/x) as X — 0; fbix) ~ g(x) as x — >■ oo; 

(iii) moment relationships fohow from 

Ef,{\X-bX-^)n = E,{Yn 

and, by R-symmetry, 

Ef,{Xn = b^+'Ef,{X-(^+^^); 

(iv) if g is decreasing, then is unimodal with mode at \/b; 

(v) if fh is unimodal, its mean and its median are both greater than its mode. 

These properties can be found in Baker [1], but only special cases of (iii) are pro- 
vided. 

Amongst the most interesting distributions presented by Baker is the root- 
reciprocal inverse Gaussian distribution (RRIG). This example, also discussed in 
Mudholkar and Wang [22], in the case of dispersion parameter A = 1, arises from 
(112.11) when g is the half-Gaussian density. The RRIG density is 

(12.2) Mx) = yie" exp |-i {x' + b^x^) | , 

& > 0. This corresponds in integral terms to (13. 3p above. Similarly, a second 
example presented by Baker |T] (Section 3.4) is the distribution based on the half- 
Subbotin distribution. This is directly linked to p.5[) . A third example, based on 
the half-i distribution, has density 

(12.3) uix) = + (X - b/x^M-^-^+^y^ 

^JvttT{v/2) 

Vjb > 0. The verification that fu.b{x) integrates to 1 can be done by using the 
integral /s in Theorem 14.11 Of course, other distributions in Baker [ij correspond 
to other integral formulae, while integral formulae in this article which have non- 
negative integrands correspond to other distributions. 

Transformation of scale densities are particularly amenable to having their skew- 
ness assessed by the asymmetry function 7(p), < p < 1, of Critchley and Jones 
[9] and provide relatively rare tractable examples thereof. Jones [14] shows that 
the asymmetry function associated with unimodal fb is 

(12.4) 7fc(p) = (^c2(p)+46-%/46) /cg(p) 

where Cg{p) = g~^(pg(^))- This shows that the Cauchy-Schlomilch transformation 
of scale always results in positively skewed distributions, with asymmetry functions 
decreasing in p, which become more skew as b decreases. For example, for the RRIG 
density (IT^ . 

(12.5) 7f,(p) = (v/2fe-logp - V2b) /V-logp. 

The extended Cauchy-Schlomilch transformation described in Section 11 also af- 
fords new transformation of scale distributions, as explored by Jones [T3]. Proba- 
bility densities of the form 

(12.6) fs{x)^g{\x~s{x)\) 
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for decreasing, onto, self-inverse s are described there. In this situation, properties 
(i)-(v) discussed above become: 

(i') fs can be defined to be S-symmetric about S-center xq since fs{x) — fs{s{x)). 
Here xq is defined by s{xq) = xq; 

(ii') /,(0) =0 and 

fs{x) ~ g{s{x)) as a; — 0; fsix) ~ g(x) as a; — )■ oo; 

(iii') moment relationships follow from 

and, by S-symmetry, 

A special case of the latter is that Ef^{s'{X)) = —1; 

(iv') if g is decreasing, then fg is unimodal with mode at xq; 

(v') if fs is unimodal and g is convex, its mean and its median are both greater 
than its mode. 

By way of example, Jones [HI briefly explored the half- Gaussian-based analogue 
of (|12.2p when s{x) is given by (|11.4p . This has probability density 

(12.7) /,(x) = y|exp|-^log2(e"--i; 

The fact that (|12.7p integrates to 1 is closely related to (|11.6p . Its asymmetry 
function has the form 

(12.8) 7.(P)- , ^— - log|cosh(av/-(logp)/2)|. 

Like ()12.5p . this asymmetry function is always positive and decreases in p; ()12.8p 
increases in a. 

The Cauchy-Schlolmilch transformation has thus motivated and triggers a new 
and promising area of work in distribution theory. 



13. Conclusions 

The Cauchy-Schlomilch transformation establishes the equality of two definite 
integrals with integrands related in a simple manner. Applying this transformation 
to a variety of well-known definite integrals yields examples that are beyond the 
current capabilities of symbolic languages. Our purpose in this paper is not only 
to present the many integrals considered here, but also to give an exposition of the 
salient points of the Cauhy-Schlomilch transformation so as to serve as motivating 
examples to explore further symbolic integration algorithms. 

Acknowledgements. The fourth author acknowledges the partial support of 
NSF-DMS 0713836. 



THE CAUCHY-SCHLOMILCH TRANSFORMATION 



19 



References 

[1] R. Baker. Probabilistic applications of the Schlomilch transformation. Comm. in Statistics- 
Theory and Methods, 37:14:2162-2176, 2008. 

[2] D. Bierens de Haan. Tables d'integrales definies. C. G. Van der Post, Amsterdam, 1st edition, 
1858. 

[3] D. Bicrcns de Haan. Expose de la theorie, des proprietes, des formules de transformation, 
et des methodes d'evaluation des integrates definies. C. G. Van der Post, Amsterdam, 1st 
edition, 1862. 

[4] D. Bierens de Haan. Nouvelles tables d'integrales definies. P. Engels, Leiden, 1st edition, 
1867. 

[5] G. Boros and V. Moll. An integral with three parameters. SIAM Rev., 40:972-980, 1998. 
[6] G. Boros and V. Moll. Irresistible Integrals. Cambridge University Press, New York, 1st 
edition, 2004. 

[7] G. Boros, V. Moll, and R. Nalam. An integral with three parameters. Part 2. Jour. Comp. 

Appl. Math., 134:113-126, 2001. 
[8] A. L. Cauchy. Sur unc formulc gcucralc relative a la transformation des integrates simples 

prises entre les limites et oo de la variable. Oeuvres completes, serie 2, Journal de I'ecole 

Polytechnique, XIX cahier, tome XIII, 516-519, 1:275-357, 1823. 
[9] F. Critchley and M. C. Jones. Asymmetry and gradient functions: density-based skewness 

and kurtosis. Scand. J. Statist., 35:415-437, 2008. 
[10] A. Erderly. Tables of Integral Transforms, volume I. McGraw-Hill, New York, 1st edition, 

1954. 

[11] A. Erderly. Tables of Integral Transforms, volume H. McGraw-Hill, New York, 1st edition, 
1954. 

[12] M. L. Glasser. A remarkable property of definite integrals. Math. Comp., 40(40) :561-563, 
1983. 

[13] I. S. Gradsliteyn and I. M. Ryzliik. Table of Integrals, Series, and Products. Edited by A. 
Jeffrey and D. Zwillinger. Academic Press, New York, 7th edition, 2007. 

[14] M. C. Jones. Distributions generated by transformation of scale using an ex- 
tended Schlomilch transformation. Submitted for publication, 2009. Available at 
http : //statistics . open . ac . uk/TechnicalReports/TechnicalReportsIntro . html 

[15] D. Kuccrovsky, E. Marchand, and R. D. Small. On the equality in distribution of the random 
variables X and g{X). Internat. J. Pure Appl. Math., 23:93-114, 2005. 

[16] D. F. Lawden. Elliptic Functions and Applications, volume 80 of Applied Mathematical Sci- 
ences. Springer- Verlag, 1989. 

[17] C. F. Lindman. Examen des nouvelles tables d'integrales definies de M. Bierens de Haan. 
P.A. Norstedt and Soner, Stockholm, 1891. 

[18] J. Liouville. Sur I'integrale 

rl jM+l/2(l _t)M-l/2 

Jo (a-|-6t-ct2)M+i 

J. Math. Pures Appl, 1:421-424, 1856. 
[19] J. Liouville. Sur I'integrale 

rl ^M+l/2(l _t)(i-l/2 

Jo (a + 6t-ct2)M+i 

Extrait d'une lettre de M. O. Schlomilch. Extrait d'une lettre de M. A. Cayley. Remarques 
de M. Liouville. J. Math. Pures Appl, 2:47-55, 1857. 
[20] A. P. Prudnikov Yu. A. Brychkov O. L Marichev. Integrals and Series. Gordon and Breach 

Science Publishers, 1992. 

[21] Z. A. Melzak. Companion to Concrete Mathematics. Mathematical techniques and various 
applications. Pure and Applied Mathematics. Wiley-Interscience Publication, 1973. 

[22] G. S. Mudholkar and H. Wang. IG-symmetry and R-symmetry: inter-relations and applica- 
tions to the inverse Gaussian theory. J. Statis. Planning Inference, 137:3655-3671, 2007. 

[23] M. Petkovsek, H. Wilf, and D. Zeilberger. A=B. A. K. Peters, Ltd., 1st edition, 1996. 

[24] O. Schlomilch. Analytische Studien. W. Engelmann, 1848. 



20 T. AMDEBERHAN, M. L. GLASSER, M. C. JONES, V. MOLL, R. POSEY, AND D. VARELA 

Department of Mathematics, Tulane University, New Orleans, LA 70118 
E-mail address: tamdeberhariSmath.tulane.edu 

Department of Mathematics and Computer Science, 
Clarkson University, Postdam, New York 13676 
E-mail address: larygaclarkson.math.edu 

Department of Mathematics and Statistics, The Open University, Milton Keynes, 
United Kingdom 

E-mail address: m.c. jonesSopen.ac.uk 

Department of Mathematics, Tulane University, New Orleans, LA 70118 
E-mail address: vhmOmath.tulane.edu 

Department of Mathematics, Baton Rouge Community College, 
Baton Rouge, LA 70806 

E-mail address: rposeySgmail.com 



DVD Academic Consulting 
E-mail address: xentricoSyahoo.com 



